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Abstract 



Recent data from heavy ion collisions at RHIC show unexpectedly large near-angle correlations that broaden longitudi- 
nally with increasing centrality. The amplitude of this ridge-like correlation rises rapidly, reaches a maximum, and then 
falls in the most central collisions. In this letter we explain how this behavior can be uniquely explained by initial-state 
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coordinate-space anisotropies converted into final-state momentum-space correlations. We propose v^/e^ part as a useful 
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1. Introduction 

The motivation for colliding heavy ions at facilities like 
the Relativistic Heavy Ion Collider (RHIC) at Brookhaven 
National Laboratory and the Large Hadron Collider (LHC) 
at CERN is to form a state of matter called the Quark 
Gluon Plasma (QGP) l]. Each of these collisions de- 
posit many TeV of energy into a region roughly the size of 
the colliding nuclei. That region is so hot and dense that 
quarks and gluons become the relevant degrees of freedom 
instead of hadrons Q. This QGP existed when the uni- 
verse was less than a microsecond old, and still extremely 
hot and dense. 

Two-particle correlations measured at RHIC reveal fea- 
tures unique to Nucleus- Nucleus collisions 0, Pair 
densities are commonly measured as a function of the az- 
imuthal angle and pseudo-rapidity difference between the 
particles (A0 and Arj respectively). While two-particle 
correlations in p+p and d+Au collisions show a narrow 
peak centered at A(j> and Arj = 0, the near-angle peak 
in Au+Au collisions broadens longitudinally and narrows 
in azimuth relative to p+p collisions. An analysis of the 
width of the peak for particles of all transverse momentum 
Pt finds the correlation extends across nearly two units of 
pseudo-rapidity When selecting higher px particles 
(e.g. pt > 2 GeV/c), the correlation extends beyond the 
acceptance of the STAR detector (Arj < 2) and perhaps 
as far as A77 = 4 as indicated by PHOBOS data This 
feature is known as the ridge. 

STAR data shows that the ridge amplitude rises rapidly 
with collision centrality Q before reaching a maximum and 
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falling off in the most central bins. The drop is often ig- 
nored. In this letter we present an explanation for the cen- 
trality dependence of the ridge amplitude related to den- 
sity inhomogeneities in the initial overlap region. We use 
measurements of the second harmonic momentum-space 
anisotropy V2 and single particle rapidity density dN/dy 
along with a Monte-Carlo Glauber model [5| for the initial 
density to predict the amplitude (Ai) of the near-side ridge 
correlation as a function of centrality for ^/snn = 62.4 
GeV, 200 GeV, and 2.76 TeV. A successful description of 
Ai is noteworthy, because it demonstrates that the corre- 
lations measured in heavy-ion collisions represent an image 
of structures such as flux-tubes @, @] in the initial overlap 
of heavy ion collisions of the order of 10 -15 m in size. 

In Glauber Monte-Carlo models @ the density dis- 
tribution in the initial collision region is assumed to be 
determined by the positions of the nucleons participat- 
ing in the interactions (participants). For a typical col- 
lision, the shape of the overlap region will be predomi- 
nantly an ellipse. The eccentricity of the ellipse can be 
quantified as e s td = fesyqifej ■ Interactions amongsts the 
systems constituents can convert that elliptic shape from 
coordinate-space into momentum-space leading to a large 
second Fourier component in the momentum-space distri- 
bution i.e. V2 = (cos 2(0 — \&rp)) > 0, where 4> is the 
azimuth angle of particles emitted from the collision and 
^rp is the reaction plane angle defined by the vector con- 
necting the centers of the two colliding nuclei. There are a 
finite number of participants in each collision so there will 
be event-to-event fluctuations in the density distributions 
that lead to fluctuations in the initial eccentricity as well as 
a tilt of the elliptic shape of the overlap region away from 
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^Rp. The plane containing the titled axis and beam axis 
is called participant plane, and the eccentricity calculated 
relative to that axis is called the participant eccentricity 
£part H • £part is a positive definite quantity so it can lead 
to non-zero even for head-on collisions. Fluctuations 
in e part can lead to fluctuations in v 2 0. It was argued 
previously that eccentricity fluctuations can give rise to 
v n fluctuations for more harmonics than just n = 2 
and that those fluctuations could therefore be the source of 
the near-side ridge [10( (especially if the v n fluctuations de- 
pend on Ar/ i.e. (v n (r])v n (r] + Arfj) = /(At?)). Subsequent 
calculations from several groups 00111 El support this. 
To test this conjecture, we calculate the centrality depen- 
dence of the near-side ridge amplitude A\ from eccentricity 
fluctuations. 



2. Eccentricity Fluctuations, Length Scales, and 
The Ridge 



We base our calculations on three premises: 1) the 
expansion of the fireball in heavy-ion collisions converts 
anisotropies from coordinate-space into momentum-space, 
2) the conversion efficiency increases with density, and 3) 
the relevant expansion plane is the participant plane. The 
participant plane can be defined for any harmonic number. 
A system with a lumpy initial energy density will lead to 
finite participant eccentricity at several harmonics 12J, |15[ ; 
eccentricity can be thought of as a harmonic decomposi- 
tion of the azimuthal dependence of the initial density. To 
illustrate how eccentricity fluctuations can lead to a ridge- 
like structure in particle correltions, using the definition 
in Ref. 0, we calculate s 2 n>pait = co S ( n ^f+(r- sin(n0)) 2 
for all harmonics. 



Fig. Q] (a) shows the n -harmonic participant eccen- 
tricity (£^p ar t) for central Au+Au collisions from our Monte- 
Carlo Glauber model. Typically the participant eccentric- 
ity is calculated based on the positions of point-like par- 
ticipants (i.e. the participant is said to exist at a precise 
x and y). One can also calculate the eccentricity from a 
more realistic model with participants smeared over some 
region. This is done by treating each participant as many 
points distributed within a disk of radius r part . The figure 
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for 7"p ar t = 0,1,2,3, and 4 fm. Increasing 
terms. The depen- 



dence of (ez 



on n is well fit with a Gaussian e 2 
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for all values of r part with the width of the Gaussian <j n 
narrowing as r part is increased. 

Here we've introduced the length scale r part causing 
the higher terms in (e^ part ) to be washed out. That effect 
is more general though and we believe it is important for 
understanding correlations and v n fluctuations [la ]. One 
can also consider what happens when particles free-stream 
for some amount of time Tf s before they interact; which 
also introduces a length scale ct{ s leading to a reduction 
of higher terms 14[. The mean- free-path (l m t p ) a particle 
travels between interactions in the fireball also affects the 
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ability of the fireball to convert higher (e^ part ) terms into 
I 3 , 17 1 . The Z m f p will prevent higher (e^ part ) terms 

from being converted into vi. The acoustic horizon may 

also play an important role |16| . In our calculations, we 
study the effect of varying r part with the understanding 
that many length-scales, like Z m f p and CTf s , will contribute 
to the final dependence of v n on n. 

is related to the two-particle correlations dN pa i r / dA(f> 
by a Fourier transform. If v„ vs. n has a Gaussian shape 
then dNp a i r / dAcf> will also have a Gaussian peak at Acp = 
0. Fig. Q] (a) shows that (e^ jPart ) follows a Gaussian so 
that eccentricity fluctuations should lead to a near-side 
Gaussian with a width depending on the length scales in 
the system. Fig [1] (b) shows the shape of the two-particle 
correlations from (e: 



n,part/ 



For n = 1, (4, P art) = since 
( e « part) ls calculated in the center-of-mass frame of the 
participants. This leads to a Gaussian centered at Atfi = 
and an apparent negative cos(A</>) term due to the sup- 
pression of the n = 1 eccentricity fluctuations. The neg- 
ative cos(A</>) is seen in the data to follow the same cen- 
trality dependence as the nearside Gaussian [3[. The fact 
provides evidence that both the nearside Gaussian and the 
awayside negative cos(A</>) terms are related to initial den- 
sity fluctuations. 

We fit the correlations in panel (b) with a Gaussian of 
width and a cos(A</>) term. In Fig. [T] (c) we plot a n 
(scaled by 1/5 for clarity) and cta<£- 0Vi(?"part) is fit with 
a power law crn^part) = 7.22r pa rt 916 and is described 
the inverse. Increasing length scales in the system should 
lead to a broadening of the azimuthal width of the near- 
side Gaussian. Having demonstrated that it is natural for 
eccentricity fluctuations to lead to a ridge like correlation 
with a width dependent on the length scales in the system, 
we now consider the amplitude of the ridge. 

3. The Ridge Amplitude 

The ridge amplitude A\ is found from data by fitting 
Apf ^/p ro f (the pair density p minus the reference pair 
density p re f scaled by ^/p rc i) vs. A(f> and Ar] The 
fit function has Ar? independent cos(A</>) and cos(2A</>) 
terms, and a near-side 2-D Gaussian with amplitude A±. 
We work with the conjecture that the 2-D Gaussian is a 



manifestation of (e„ part ) and calculate the centrality de- 
pendence of A\. Our result for A\ will be related to so 
we need to know the transfer function or conversion effi- 



ciency c„ [l6| of (e 



n.part, 



into v^- We expect c„ to depend 
on particle density. In Fig. [5] we show the two-particle 
cumulant scaled by eccentricity f2 |2|/ g2. p art vs. density 
{l/S)dN/dy from RHIC and LHC [ll [20j] . The funciton 
from Ref. [2l| fits the full range of RHIC and LHC data 
reasonably well providing our estimate of C2 . 

We can estimate A\ independently from the fluctua- 
tions of several different harmonics. Since a cos(2A<^) term 
proportional to £g td is used in the fit to Ap/^/p ro f [22j we 
base our estimate of n — 2 component of A\ on the differ- 
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We can also calculate the 
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Figure 1: a) (e^ part ) for central Au+Au collisions from a Glauber Monte-Carlo where participants are treated as point-like or smeared over 
a region of size r par t. The curves show Gaussian fit functions, b) the A</> dependence of two-particle correlations if (v^) <x (e^ part }- The 
curves show a Gaussian+cos(A</>) fitting function. Bottom panel: The Gaussian widths extracted from the fits in panel a) and b) with a 
power-law fit to describe the r pa rt dependence. 
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Figure 2: The ratio of i>2{2} over £2, part from a Glauber Monte- 
Carlo vs. (l/S)dN/dy for = 0.2 and 2.76 TeV. The fit function 
^{2}/£ 2 ,part = A/(l + (1/s) B dN/dy ) is adapted from Ref. gj]. The 
fitting parameters are A = 0.412 and B = 15.175. 



n = 3 component from e| part alone but this requires an 
assumption about C3 which is not yet measured. We need 
to convert our prediction for correlations from a particu- 
lar harmonic into a Gaussian amplitude A\. For n = 2 we 
find 
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2^ 



Ap 



cos(2A</>)dA</> = O.llAi 



(1) 



with the n = 3 component similary giving 0.039^4i. The 
azimuthal width cta^ of the near-side Gaussian is weakly 
dependent on centrality so we use a typical value of cta^ = 

'p ie f we need to include the 



0.65. To relate v 2 n to Ap/ V 
particle density p = ±^§- 



Combining (e| part ) 



( £ std) w ith the conversion efficiency C2, particle density po, 



and factor of 0.11 we find A\ poC2\s 
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2. part 
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)/0.11. 



We take po from data, c 2 from Fig. [2j and (e^part) an d 
£ s td from our Monte-Carlo Glauber model. The n = 3 
based prediction does not need e s td but does require the 



conversion efficiency C3. We assume that C3 has the same 
density dependence as c-i but that various length scales in 
the system will suppress C3 relative to C2. C3 = C2/2.6 gives 
good agreement between the n = 2 and n = 3 estimates 
and is consistent with the AMPT results in Ref. [l2| . From 
n = 3 we find A\ ~ po c 3£|,part/0-039 where C3 ss C2/2.6. 

Fig. [3] shows our estimate of the ridge amplitude A± 
based on (e^ paxt ) vs. centrality parameter v = 2A r bi n /^ r part 
for Au+Au collisions at 62.4 and 200 GeV and for Pb+Pb 
collisions at 2.76 TeV. iVbin and A part are the number of 
binary nucleon-nucleon collisions and the number of par- 
ticipants. More central or higher energy collisions yield 
larger values of v. The open (closed) symbols show the 
n = 3 (n = 2) estimates. 
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Figure 3: An estimate of the contribution of eccentricity fluctuations 
to the near-side Gaussian peak amplitude A\ vs. centrality measure 
v = 2N Mn /N pa rt for 62.4 Au+Au, 200 GeV Au+Au, and 2.76 TeV 
Pb+Pb collisions. The closed symbols are the amplitudes expected 
from £j par t — £ s t d an d t ne closed symbols are from e| paxt - 
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4. The Rise and Fall 

We find that our estimate of Ai agrees with what has 
been observed at 200 and 62.4 GeV. Our A\, like the data, 
starts at a small value and rises much faster than expecta- 
tions from a linear superposition of independent p+p colli- 
sions (shown for 200 GeV in the Fig. [3]) which assumes that 
correlations grow as Nbi n /dN/dy. The rise continues until 
Ai reaches a maximum of 0.7 for 200 GeV near v = 5, then 
falls again. The 2.76 TeV calculations show a similar trend 
as the 200 and 62.4 GeV calculations but the amplitude is 
expected to be much larger. This provides a testable pre- 
diction for the LHC experiments. Our picture provides a 
natural explanation for the rise and fall related to the ini- 
tial overlap geometry and its fluctuations. For the n = 3 
term for example, (£l iPart ) falls with -/V part since the larger 
iVp ar t leads to smaller fluctuations. But iVpart(£§ par t) nrs ^ 
rises then falls (see Fig. 2]). This rise and fall is due to the 
almond shaped geometry of the overlap region which leads 
to non-statistical fluctuations in e\ part for more than just 
the n = 2 harmonic [15j . Fig. [j shows N part e\ for 
n = 3, 5, 7, 9, and 23. For harmonics close to n = 2, the 
almond shape of the overlap zone causes large deviations 
from a trivial £„ part cx 1/N part behavior. The higher har- 
monics approach this statistical expectation with n = 23, 
nearly reaching a 1/N part behavior. 
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Figure 4: The centrality dependence of odd terms for N par t^n part' 
The rise and fall of the ridge can be traced to the non-statistical 
centrality dependence of e\ part caused by the intrinsic almond shape 
of the nuclear overlap region. 

Both c n and po are increasing with centrality but the 
product of /9oc n (£n P art) rises until very central collisions 
and then falls. The drop fall in central is due to the sub- 
sidence of the effect of the almond shape of the overlap 
geometry. The observation that the rise and fall exists in 
the near-side ridge amplitude shows that the ridge is dom- 
inated by initial geometry fluctuations; even exhibiting a 
dependence on the coupling of various harmonics to the in- 
trinsic almond shape of the overlap region. We know of no 



other plausible scenarios to explain the rise and fall of the 
ridge other than this explanation. This rise and fall is to 
our knowledge a unique signature of density fluctuations. 
This implies that momentum-space correlations are sensi- 
tive to initial density fluctuations and that correlations in 
heavy-ion collisions provide an image of the density distri- 
butions in the initial overlap zone. The system created in 
these collisions acts as a femtoscope, revealing structures 
with length scales on the order of a fm. 

5. Discussion and Conclusion 

Our explanation for the centrality dependence of A\ 
provides a natural explanation for the rise and fall of the 
ridge. Our estimates of the amplitude agree with RHIC 
data suggesting that the conversion of geometry fluctua- 
tions in the initial overlap region into momentum space 
causes the near-side ridge structure. We expect the con- 
version efficiency to drop with n since effects like initial- 
state free-streaming and mean-free-path will wash out the 
higher harmonic terms. Extracting the conversion effi- 
ciency c n = - ""^| 2 } 2 as a function of n [l6j], centrality, 
and particle kinematics will provide information on those 
effects. This only relies on measuring the two-particle cu- 
mulant v n {2} 2 which is a rather simple experimental mea- 
surement and comparing it to the initial eccentricities from 
various models of the initial density. It will be particularly 
interesting to determine c„ as a function of Arj to under- 
stand how de-coherence affects manifest in the longitudinal 
direction. 

We presented the participant eccentricity vs. harmonic 
when the participants are treated as point-like or smeared 
over a radius r part . The larger values of r part wash out 
the higher harmonic eccentricities. We argued that, sim- 
ilarly, a large mean-free-path or acoustic horizon should 
also wash out higher harmonics of v n . Such an effect 
could lead to a Gaussian peak in two-particle correlations 
at Acf> = similar to that seen in the data. We calcu- 
lated the contribution to the near-side Gaussian peak that 
we expect from initial density fluctuations. Following sim- 
ple premises, we find that the near-side peak from density 
fluctuations should rise rapidly, reach a maximum just be- 
fore the most central events, then fall. Our estimate of the 
magnitude is in agreement within our uncertainties with 
the available data and the shape matches the data. This 
is the only calculation we know of to correctly describe the 
rise and fall of the ridge amplitude. Our calculation shows 
that the rise and fall is related to the interaction of fluc- 
tuations with the shape of the initial overlap geometry. 
We conclude therefore that density fluctuations are the 
dominant source for the ridge-like correlations. The longi- 
tudinal width of the ridge remains an interesting topic to 
investigate. Finally we use the same framework to predict 
the ridge amplitude for 2.76 TeV Pb+Pb collisions. 

The authors thank Sergei Voloshin, Sean Gavin and 
Joern Putschke for their helpful comments. 
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